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ABSTRACT 

In QCD with two flavors of massless quarks, tlie cliiral pliase transition is 
plausibly in the same universality class as the classical four component Heisen- 
berg antiferromagnet. Therefore, renormalization group techniques developed in 
the study of phase transitions can be applied to calculate the critical exponents 
which characterize the scaling behaviour of universal quantities near the critical 
point. This approach to the QCD phase transition has implications both for lattice 
gauge theory and for heavy ion collisions. Future lattice simulations with longer 
correlation lengths will be able to measure the various exponents and the equa- 
tion of state for the order parameter as a function of temperature and quark mass 
which we describe. In a heavy ion collision, the consequence of a long correlation 
length would be large fluctuations in the number ratio of neutral to charged pions. 
Unfortunately, we show that this phenomenon will not occur if the plasma stays 
close to equilibrium as it cools. If the transition is far out of equilibrium and can 
be modelled as a quench, it is possible that large volumes of the plasma with the 
pion field correlated will develop, with dramatic phenomenological consequences. 



1. Introduction 

The QCD phase transition is of interest from several different points of view. 
First, there can be no doubt that it occurred in the early universe. Second, it 
is reasonable to hope that in a heavy ion collision of sufficiently high energy, a 
small region of the high temperature phase is created which then cools through the 
phase transition. Third, lattice gauge theory is well suited to calculating the equi- 
librium properties of QCD at high temperatures. From all these perspectives, it is 
important to learn as much as can be learned analytically about the phase transi- 
tion, relying as much as possible only on fundamental symmetries and universality 
arguments and as little as possible on specific assumptions and models. 

In a previous paper [1], one of us (F. W.) emphasized that in the chiral limit 
where there are two species of quarks with zero current algebra mass, the order pa- 
rameter for the chiral phase transition has the same symmetry as the magnetization 
of a four component Heisenberg magnet, which has a second order phase transi- 
tion. In this case, and indeed for any number of quark species except zero, there 
is no order parameter for a confinement /deconfinement phase transition. Thus in 
the niu = m^ = limit the universal characteristics of the QCD phase transition 
{i.e. those characteristics determined by the modes which develop long correlation 
lengths at the phase transition) are the same as those of the A^ = 4 Heisenberg 
magnet. 

In this paper we further explore the consequences of this approach to the QCD 
phase transition. In the following section, we review the scenario described in [1], 
and establish a dictionary between QCD and the magnetic system. In order to 
make the present paper self-contained, we also include in this section many of the 
results from [1]. In section 3, we discuss the behaviour of the pion and sigma 
masses at the transition. In section 4, we discuss the p and Ai mesons. In sec- 
tion 5, we consider how the strange quark affects the phase transition. In section 
6, we go beyond the static critical phenomena of the earlier sections and discuss 
the dynamics of the appropriate universality class. In section 7, we discuss the 



implications of all this for cosmology and lattice gauge theory, and for heavy ion 
collisions under the assumption that the plasma remains close to thermal equilib- 
rium through the phase transition. In section 8, we consider what phenomena we 
can expect in heavy ion collisions if the system gets far out of thermal equilibrium 
and can be modelled as a quench. Finally, in the last section we summarize and 
conclude. 

2. QCD and the 0(4) Magnet 

As discussed in [1], the physics of the QCD phase transition is qualitatively 
different in the cases of zero, one, two, or three or more flavors of quarks. In this 
section we consider QCD with two species of quarks. (An analysis similar to the 
one which follows leads to the conclusion that for three or more flavors of massless 
quarks, the chiral phase transition is first order. See [1] for details.) If there are 
two fiavors of massless quarks, the lagrangian is symmetric under global chiral 
transformations in the group SU{2)l x SU{2)^ x f/(l)i+i? of independent special 
unitary transformations of the left and right handed quark fields, and a vector 
f/(l) transformation which corresponds to baryon number symmetry. (The axial 
f/(l) which would make the symmetry group into U{2) x f/(2) is a symmetry of the 
classical theory, but not of the quantum theory [2].) This chiral symmetry breaks 
spontaneously down to SU(2)i^^ x f/(l)L+i? at low temperatures; and is restored 
at sufficiently high temperatures. The order parameter for this phase transition is 
the expectation value of the quark bilinear 

M) ^ (qlqRj) (2.1) 

which breaks the symmetry when it acquires a non-zero value below some critical 
temperature Tc. 

In order to describe a second-order transition quantitatively, we must find a 
tractable model in the same universality class. For the chiral order parameter (2.1) 



the relevant symmetries are independent unitary transformations of the left- and 
right-handed quark fields, under which 

M -^ U^MV . (2.2) 

These transformations generate an U{2)i x U{2)^ symmetry, which is not quite 
what is needed, since it includes the axial baryon number symmetry which is not 
present in QCD. This problem is solved [1] by restricting Ai to unitary matrices 
with positive determinant, instead of general complex matrices. Matrices Ai in 
this restricted class remain in this restricted class under the transformation (2.2) 
only if U and V have equal phases. Hence, the axial U{1) has indeed been removed. 
The 2x2 matrices Ai can conveniently be parametrized using four real parameters 
(o", n) and the Pauli matrices as 

M = a + irr -r . (2.3) 

In fact the order parameter can be written as a four-component vector (p = {a, tx) 
and the transformations (2.2) are simply 0(4) rotations in internal space. Hence, 
the order parameter appropriate for the chiral phase transition in QCD with two 
flavors of massless quarks has the symmetries of the standard 0(4) invariant A^ = 4 
Heisenberg magnet. For smaller number of components, this sort of model is a 
much-studied model for the critical behavior of magnets, with the order parameter 
representing the magnetization of a ferromagnet or the staggered magnetization of 
an antiferro magnet. 

If the phase transition is second order, then it will correspond to an infrared 
fixed point of the renormalization group. We wish to describe those aspects of 
the critical behaviour which are universal, that is, those aspects which are deter- 
mined by the scaling behaviour of operators near the infrared fixed point of the 
renormalization group. Hence, it is sufficient to retain those degrees of freedom 
which develop large correlation lengths at the critical point. These are just long 



wavelength fluctuations of the order parameter, which is small in magnitude near 
the critical point and therefore fluctuates at little cost in energy. Thus the most 
plausible starting point for analyzing the critical behavior of a second-order phase 
transition in QCD is the Landau- Ginzburg free energy 

d^x{^ d'rdi(l)a + y 0>« + ^(0>a)2 } . (2.4) 

Here /i^ is the temperature-dependent renormalized (mass)^, which is negative 
below and positive above the critical point, while A is the strength of the quartic 
couplings and is supposed to be smooth at the transition. We neglect terms with 
higher powers of since |0| is small near the transition. The symmetry breaking 
pattern we want is A1 oc 1 (equivalently, (o") 7^ 0; (vf) = 0) below the transition 
which is indeed what we find at the minimum of the potential for positive A. This 
model has been studied in depth for arbitrary A^ and spatial dimension d, and the 
existence of an infrared stable fixed point of the renormalization group has been 
established [3]. Hence, it is a model for a second order QCD chiral phase transition 
for two massless quarks. 

When the free energy (2.4) is written in terms of a and vf it looks much like the 
original model of Gell-Mann and Levy [4] with two changes: there are no nucleon 
fields and only three (spatial) dimensions. These two changes refiect an impor- 
tant distinction [5]. We are only proposing (2.4) as appropriate near the second 
order phase transition point. This is because it is only there that we can appeal to 
universality - the long-wavelength behaviour of the a and vf fields is determined 
by the infrared fixed point of the renormalization group, and microscopic consid- 
erations are irrelevant to it. In Euclidean field theory at finite temperature, the 
integral over u of zero temperature field theory is replaced by a sum over Matsub- 
ara frequencies Un given by 2miT for bosons and (2n+ l)7rT for fermions with n an 
integer. Hence, one is left with a Euclidean theory in three spatial dimensions with 
massless fields from the n = terms in the boson sums and massive fields from the 
rest of the boson sums and the fermion sums. Hence, to discuss the massless modes 



of interest at the critical point, (2.4) is sufficient. This means that we do not need 
to worry whether to introduce nucleon fields as in [4], or constituent quark fields 
as, for example, in [6]. 

We have motivated a very definite hypothesis for the nature of the phase tran- 
sition for QCD with two species of massless quarks, namely that it is in the uni- 
versality class of the N = 4 Heisenberg magnet. This hypothesis has numerous 
consequences, which are the subject of the rest of this paper. To keep the discus- 
sion self-contained, in the remainder of this section we review the predictions for 
the static critical exponents described in [1]. 

2.1. Critical exponents 

First, we define the reduced temperature t = (T — Tc)/Tc. The exponents 
a, (3, 7, 7], and u describe the singular behaviour of the theory with strictly zero 
quark masses as t — > 0. For the specific heat one finds ce[6 

C{T) ~ Itp" + less singular. (2.5) 

The behaviour of the order parameter defines (3. 

- \tf for t < . (2.6) 



1] and u describe the behaviour of the correlation length ^ where 

Gaf3{x) = (0(a;)a0(O)/3) - (0a)(0/3) ^ 5^/3— exp(-|x|/0 at large distances. 

(2.7) 
A is independent of |a;|, but may depend on t. The correlation length exponent z/ 
is defined by 

e~itr"- (2.8) 

Above Tc, where the correlation lengths are equal in the sigma and pion channels. 



the susceptibility exponent 7 is defined by 

c^\G,^(x)~t-^ (2.9) 

We will discuss the behaviour of the susceptibility below the transition in the 
following section. The exponent rj is defined through the behaviour of the Fourier 
transform of the correlation function: 

Gapik ^ 0) ~ A;-2+^ . (2.10) 

The last exponent, S is related to the behaviour of the system in a small 
magnetic field H which explicitly breaks the 0(4) symmetry. Let us first show 
that in a QCD context, H is proportional to a common quark mass rriu = rrid = 
niq. This common mass term may be represented by a 2 x 2 matrix V given by 
iTLq times the identity matrix. We are now allowed to construct the free energy 
from invariants involving both V and M.. The lowest dimension term linear in 
T> is just iiJ^/vV = niqa, which in magnet language is simply the coupling of the 
magnetization to an external field H oc rriq. In the presence of an external field, 
the order parameter is not zero at T^. In fact, 

m)it = 0,H^O)r^H'/' . (2.11) 

The six critical exponents defined above are related by four scaling relations 
[3]. These are 

a = 2 — dv 

i5 = ^(rf-2 + r^) 

7 = {2-r])u ^^-^^^ 

d + 2-ri 
d-2 + r] ■ 

We therefore need values for 77 and z/ for the four component magnet in li = 3. 
These were obtained in the remarkable work of Baker, Meiron and Nickel [7], who 



carried the perturbation theory to seven-loop order, and used information about 
the behaviour of asymptotically large orders, and conformal mapping and Pade 
approximant techniques to obtain 

7] = .03 ±.01 

(2.13) 
u = .73 ±.02 . 

Using (2.12) , the remaining exponents are a = —0.19 ± .06, (3 = 0.38 ± .01, 
7 = 1.44 ± .04 and 6 = 4.82 ± .05. Since a is negative there is a cusp in the specific 
heat at T^, rather than a divergence. 

Very different critical exponents are proposed in [8]. These authors thus implic- 
itly claim that a hitherto unknown fixed point theory with the symmetries of the 
N = 4 isotropic Heisenberg model exists, and governs the QCD phase transition. 

3. The Equation of State and the Pion and Sigma Masses 

The expressions which define (3, 7 and 6 are actually special cases of a more 
general relationship between the magnetization and the magnetic field called the 
critical equation of state. The equation of state has been calculated to order e^ 
by Brezin, Wallace and Wilson [9]. Their result is reproduced in the appendix. In 
this section, we will use the equation of state to determine the behaviour of the 
masses of the pion and sigma masses near the critical point. 

First, we must define what we mean by the "mass" of the pion and sigma. We 
could choose either to define the mass as an inverse correlation length or as an 
inverse susceptibility. We choose the latter, which is conventional in the condensed 
matter literature. Specifically, we define 

m^2 ^ I d^xGoo (3.1) 

and 

m~Hij = f d^xGij (3.2) 

where (po = a and (pi = TTi,i = 1,2,3. At any given t and H, (2.7) implies that 



whether one defines the mass as the inverse correlation length or as the inverse 
susceptibility is academic. However, since A in (2.7) depends on t and H, the two 
different choices lead to different scaling behaviours for masses as functions of t 
and H. We shall see that with the conventional choice (3.1) and (3.2) , the masses 
can be extracted conveniently from the equation of state. It is worth noting that 
the masses we have defined are related only to the behaviour of spatial correlation 
functions in the static (equilibrium) theory. They carry no dynamical information. 
Also, we will only be able to make universal statements about how the masses scale 
at the transition. Normalizing the magnitudes of the masses {i.e. relating them 
to the zero temperature masses) will require using some specific model, and hence 
will not be universal. 

The equation of state gives the magnetization as a function of t and H. For 
the rest of this paper, we will write the order parameter as M, for magnetization, 
keeping in mind that M = (a) = (|0|). In order to define the equation of state, 
we first define a shifted field a = a — {a) = a — M. Then the equation of state is 
simply the relation 

{a) = 0. (3.3) 

This relation has been expanded to order e^ by Brezin, Wallace and Wilson [9]. 
The result can be expressed conveniently in terms of the variables y = H/M and 

X = t/M^II^ as 

y = fix) (3.4) 

where the function f(x) was calculated to order e^ in [9], and is given in the 
appendix. The units in which H and M are measured are chosen so that /(O) = 1 
corresponds to t = and /(—I) = corresponds to the t < 0, H = coexistence 
curve. Knowing f{x), we can calculate the value of the order parameter M for a 
given H and t using (3.4) . The behaviour of the order parameter is illustrated 
in Figure 1. This figure and the other ones in this section should be viewed as 
illustrations of qualitative behaviour rather than quantitative predictions because 
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they are based on setting e = 1 in the 0(e) expression for f{x). The values 
for the critical exponents themselves which we quoted in the previous section are 
quantitative predictions, complete with error estimates, because they are based on 
the much more elaborate analysis of Baker et al. [7] 

From the equation of state, we can deduce the behaviour of m^r and nia at 
non-zero (but small) t and H. The masses are given by 

o dH , , 

and 

'nl = |. (3^6) 

The first relation follows directly from the definition (3.1) , and the second follows 
from (3.2) and from assuming that M || H, so that a small change 6H ± H gives 
a small change 6M ± M with 6M/6H = M/H. Using the equation of state, we 
can rewrite (3.5) and (3.6) as 

ml = M^-'fix) (3.7) 

and 



ml = M^-i ( 6 fix) - -J'{x) ) . (3.^ 



Hence, as suggested in [8], 5 can be determined by measuring the ratio mil /ml at 
t = 0. In general, from f{x) we can find the pion and sigma masses for any t and 
H. 

There are two interesting limits which we will consider explicitly. First, for 
t > and H ^ which corresponds to a; ^ oo, we should find the full 0(4) 
symmetry, and hence should find that the pion and sigma masses are identical. 
For X -^ oo, the function f{x) from the appendix behaves as f{x) = cx'~^ . Here, 



11 



the constant c and the exponent 7 are given to 0(e) and 0(e) respectively in (A9) 
and (AlO) . Applying (3.7) and (3.8) , we find that 

m^ = m^ = cf^ for a; — !► 00 , (3.9) 

consistent with the symmetry. 

We can also consider the limiting case of approaching the coexistence curve. 
This means taking t < and H ^ 0, which implies x — ^ — 1. In this limit, M tends 
to a nonzero constant, and so from (3.6) , we obtain m^ oc H, a, familiar result for 
Goldstone bosons. The behaviour of the pion mass is illustrated in Figure 2a. The 
result (3.6) may look peculiar to a particle physicist who is more familiar with the 
zero temperature result 

2 _ 2mq{qq) 
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J -K 
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(3.10) 



Before considering the sigma mass, we therefore pause here to explain how (3.10) 
and (3.6) are related. We have seen that niq ~ H and that the order parameter 
(gg) ~ (o") ~ M . At zero temperature, f-,^ is defined in terms of the axial current 
by the relation 

(OM-(0)|7r^(g))=^Ag^(5"^ (3.11) 

In the zero temperature linear sigma model, the axial current is given by 

A^(x) = a(x)(9^7r°(x) - 7r"(x)(9^(T(a:) , (3.12) 

which means that /tt defined in (3.10) is simply 

/7r = (0|a|0) . (3.13) 

This result suggests that we make the identification /jr ~ M, which does indeed 
make (3.6) and (3.10) equivalent. However, it is important to remember that using 
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the linear sigma model at zero temperature can not be justified by a universality 
argument in the way that using it near T = T^ can. Hence the argument of this 
paragraph is not a derivation of (3.6) from the zero temperature result (3.10) . 
(3.6) is valid near T = T^ while (3.10) is valid at T = 0. Also, m^r in (3.10) is 
a mass in a 3 + 1 dimensional Lorentz invariant theory, while mjr in (3.6) is an 
inverse susceptibility in a 3 dimensional theory. We have simply shown that a 
reader familiar with one expression should not be surprised by the other. 

The behaviour of the sigma mass at the coexistence curve is trickier to obtain 
than that of the pion mass. First, we note that in mean field theory (e = 0) the 
equation of state is simply y = f{x) = 1 + x, and rria is easily evaluated using (3.8) 
. For H ^ at fixed t < the result is 



Hence, in mean field theory m^ decreases with H to a, non-zero value at H = 0. 
However, for d < A when fiuctuations are important, the result is quite different. 
In words, fluctuations of the massless pions produce new infrared singularities in 
the longitudinal susceptibility, or, equivalently, make the sigma massless. Now, let 
us see how this result can be obtained from the equation of state [10,11]. In the 
limit H ^ 0, f{x) ~ H while f'{x), we will see, tends to zero more slowly. Hence, 
the second term in (3.8) is dominant and gives 

f3ml 



ikf^-i 



f\x) for X ^ -1 . (3.15) 



The difficulty is that from the expression (All) for f{x) valid for a; — ;► — 1, we notice 
that f'{x) contains divergent terms like e log(a; + 1), e^ log (a; + 1) and e^ log(a; + 1). 
These terms do not exponentiate to f'{x) ~ (a; + 1)^, contrary to the claims of [8]. 
After some algebra [11], one flnds the result 

Both the terms on the right side of (3.16) must be kept because they differ in their 
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exponents only by order e. Also for this reason, the constants ci and C2 given in 
the appendix are known only to order e even though f{x) is known to order e^. 

Qualitatively, as H is lowered at fixed t < 0, at first the ci term dominates 
and m^ appears to be decreasing toward a non-zero value at /J = as in the mean 
field result. Then, the C2 term takes over and one finds that in fact the sigma mass 
goes to zero like m^ oc H'^' . The behaviour of the sigma meson mass is illustrated 
in Figures 2b and 3. In future lattice simulations, as rriq is lowered toward zero, 
this behaviour should be observed. This result is an example of the power of the 
renormalization group techniques in obtaining universal results. If we had chosen 
a specific model, say that of Gocksch [6], or the Nambu-Jona-Lasinio model of 
Hatsuda and Kunihiro [12], we would have been able to calculate non-universal 
quantities far from T^, but would basically have been limited to using mean field 
theory, as those authors do. Then, we would have reached the incorrect conclusion 
that the sigma has a non-zero mass in the chiral limit below Tc- Here, by restricting 
ourselves to calculating universal quantities, we are limited to the region near the 
critical point, but we can be confident in our results regardless of which specific 
model is correct and can include the effect of fiuctuations. 

4. p and al Correlation Functions 

To this point, we have discussed the correlation functions in the pion and sigma 
channels only. It is certainly possible to construct other spatial correlation func- 
tions. The next-simplest are those associated with the p (Lorentz vector, isospin 
vector) and Ai (Lorentz vector, isospin axial vector) mesons. In Gocksch's model 
[6], these have correlation lengths given simply by 2mQ{T) where rriQ = iiT + g{a) 
is the mass of the constituent quarks in his model. Since these constituent quarks 
are fermions, they have a Matsubara mass of vrT. g is a, coupling constant. In any 
model, there is bound to be a model dependent contribution to these correlation 
lengths which is smooth at Tc- However, the pions and sigma also contribute to 
the p and Ai correlation lengths since operators with the appropriate symmetries 
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can be constructed from the pion and sigma operators. In particular, 

P^ = ea/377r/3VV^ (4.1) 

and 

(Ai)^ = aVV,-7r,VV . (4.2) 

Hence, the pions and sigma can make a universal non-analytic contribution to the 
p and Ai correlation lengths at the critical point. In this section, we calculate this 
contribution for t > and H = 0. 

Although we will try to define all the quantities we use, for those interested in 
further details we note that we are following the notation and conventions of Amit 
[3]. In the region t > and H = the symmetry is unbroken and hence there is 
no distinction between the a and tt operators which we will call 0" and the p and 
Ai operators which we combine into Of = e°'^^ (jf^V i(j) . We are interested in 
the scaling behaviour of the correlation function {0{x)0{Q)) = r*- ' ' -*. In general, 
by r*-"^'"''^-' we mean the m-point vertex function with n insertions of and p 
insertions of O. The scaling behaviour of vertex functions involving a composite 
operator like O is determined by rj, u, and the anomalous dimension of the operator. 
Here we are fortunate because the operators Oi are the "conserved" currents of 
the chiral symmetry, where in the 3 dimensional theory this means V^Oi = 0. The 
Ward identities arising from chiral symmetry imply that the O are not subject to 
renormalization [13], or, equivalently, that their anomalous dimensions are zero. 

For H = 0, the vertex function r'^''^'^-* is a function of the external momentum 
k, the reduced temperature t, the renormalization point k, and the (p^ coupling 
constant A. At the fixed point of the renormalization group, A is a constant and 
the vertex function satisfies the renormalization group equation [3] 

;|:-(^^)4}^'"'°'^*(*-*-^'-«' <^-^' 

The right hand side of the equation is non-zero because although T^'^'^''^' is not 
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subject to multiplicative renormalization, it is additively renormalized. However, 
B does not depend on either k or t [3], and therefore does not contribute to the 
no n- analytic behaviour at the critical point. Therefore, in what follows we will 
drop B. As a consequence of (4.3) , the vertex function has the form 

T^'''''\k,t,K)=F(^K'/\ k) , (4.4) 

where 9 = - — 2. By dimensional analysis, we know that 

This means that (4.4) becomes 

r<°*^' = ''"f(f(^)"".j). («) 

To this point, b has been arbitrary. Therefore, we can choose it strategically. With 
the choice b = K{t/ k^Y ^ we find 

r(o,o,2) ^ iKd-2)^(^^-.) (4 7) 

where the function g{x) satisfies g{Q) = constant. If, on the other hand, we choose 
b = k, the result is 

p(0,0,2) _ ^d-2^(^^-^) (48) 

where the function g{x) tends to a constant for a; -^ oo. (4.7) , which gives 
the behaviour of r*-^'*^'^-* as a function of t for k = 0, and (4.8) , which gives 
the behaviour for t = as a function of k, describe the non-analytic part of the 
correlation function. As we mentioned earlier, there will also be a model dependent 
but analytic mass term for the p and Ai. Except near the critical point, this 
smooth term is presumably larger than the non-analytic term whose effects we 
have calculated. In order to observe (4.7) and (4.8) , future lattice simulations will 
have to get close enough to the critical point that the non-analytic term dominates 
the analytic term. 

16 



5. The Influence of the Strange Quark 

To this point in tliis paper, we have described a world with two massless quarks, 
and hence we have imphcitly been taking the strange quark mass to be infinite. If 
the strange quark is massless, then Pisarski and Wilczek showed [1] that the chiral 
phase transition is first order. Hence, as the strange quark mass is reduced from 
infinite to zero, at some point the phase transition must change from second order 
to first order. This point is called a tricritical point. There is numerical evidence 
[14] that when the strange quark has its physical mass, the transition is second 
order. Hence, we devote most of this paper to analyzing the second order phase 
transition. However, in a lattice simulation, the strange quark mass could be tuned 
to just the right value to reach the tricritical point. In this section, we discuss the 
critical exponents that would be observed in such a simulation. 

Let us consider the effect of adding a massive but not infinitely massive strange 
quark to the two flavor theory. This will not introduce any new fields which 
become massless at Tc, and so the arguments leading to the free energy (2.4) 
are still valid. The only effect of the strange quark, then, is to renormalize the 
couplings. Renormalizing /x^ simply shifts T^, as does renormalizing A unless A 
becomes negative. In that case, one can no longer truncate the Landau-Ginzburg 
free energy at fourth order. After adding a sixth order term, the free energy 
becomes 

d'x{lwf + ^02 + A^^2)2 ^ -(^2)3 _ ^^1 . (5.1) 

While for positive A, 0^ increases continuously from zero as /i^ goes through zero, 
for negative A, 0^ jumps discontinuously from zero to \\\/{2k) when /i^ goes 
through \'^/(4:k). Hence, the phase transition has become first order. Thus at 
the value of rris where A = 0, the phase transition changes continuously from 
second order to first order. 

The singularities of thermodynamic functions near tricritical points, like the 
singularities near ordinary critical points, are universal. Hence, it is natural to 
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propose [1] that QCD with two niassless flavors of quarks and with T near Tc and 
rris near its tricritical value is in the universaUty class of the 0^ Landau-Ginzburg 
model (5.1) . This model has been studied extensively [15]. Because the inter- 
action is strictly renormalizable in three dimensions, this model is much simpler 
to analyze than the 0^ model of the ordinary critical point. No e expansion is 
necessary, and the critical exponents all take their mean field values. There are 
calculable logarithmic corrections to the scaling behaviour of thermodynamic func- 
tions [15], but we will limit ourselves here to determining the mean field tricritical 
exponents. 

In mean field theory, the correlation function in momentum space is simply 
Gapik) = 5ap{k'^ + /i^)~^. Since /x^ ~ t, this gives the exponents r^ = 0, 7 = 1 and 
V = 1/2. To calculate a and [3, we minimize F for H = X = V0 = 0, and find 
a = 1/2 and (3 = 1/4. To calculate 6, we minimize F for t = X = V0 = and find 
6 = 5. 

The result for the specific heat exponent a is particularly interesting, since it 
means that the specific heat diverges at the tricritical point, unlike at the ordinary 
critical point. This means that whereas for rris large enough that the transition is 
second order the specific heat C{T) has a cusp but is finite at T = Tc, as rris is 
lowered to the tricritical value C{Tc) should increase since at the tricritical point 
it diverges. This behaviour should be seen in future lattice simulations. 

Finally, at a tricritical point there is one more relevant operator than at a 
critical point, since two physical quantities (t and rris) must be tuned to reach a 
tricritical point. Hence, a new exponent 0t, the crossover exponent, is required. 
For A 7^ 0, tricritical behaviour will be seen only for \t\ > t*, while for \t\ < t*, 
either ordinary critical behaviour or first order behaviour (depending on the sign 
of A) results, t* depends on A according to 

t* ~ A^/*^* (5.2) 

The mean field value of (pt is obtained by minimizing the free energy F for H = 
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V0 = 0, and is 0t = 1/2. These mean field tricritical exponents, a = 1/2, j3 = 1/4, 
'-f = 1, 6 = 5, 1] = 0, u = 1/2, and (pt = 1/2 would describe the real world if rris 
were smaller than it is, and will describe future lattice simulations with rris chosen 
appropriately. 

6. Dynamic Critical Phenomena 

To this point, we have discussed the static critical phenomena appropriate 
to the equilibrium properties of the QCD plasma near its phase transition. To 
discuss dynamical phenomena, we need equations of motion. In zero temperature 
scalar field theory, Lorentz invariance requires that these equations have no first 
order time derivatives. In the finite temperature theory appropriate for discussing 
critical behaviour, however, there is no Lorentz invariance, and hence we should 
expect first order time derivatives. In analogy with the notion of static universality, 
one finds that the dynamics of the long wavelength modes groups theories into 
dynamic universality classes containing theories described by the same equations 
of motion. To specify the dynamic universality class, one needs to specify more 
than the dimension of space and the number of components of the order parameter. 
Hence, there are usually several different dynamic universality classes which are 
all in the same static universality class. In particular, it is necessary to specify 
whether or not the order parameter is conserved, and which other quantities are 
conserved. Thus a ferromagnet will have very different dynamical behaviour from 
an antiferromagnet with the same number of components, even though the static 
universality class is the same. The ferromagnet will have a much more difficult 
time thermalizing long-wavelength fiuctuations, since in the A; ^ limit they are 
rigorously stable. 

From this point of view, two fiavor QCD behaves as an antiferromagnet. Its 
order parameter, the expectation value of a scalar or pseudoscalar quark bilinear, 
is not a conserved quantity. (It is quantities of the form q'~fQq, not qq, that are con- 
served). This means that a model similar to model G of Halperin and Hohenberg 
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[16], is appropriate. Model G is formulated for a three component order parameter. 
We show below how to formulate it for the A^ = 4 order parameter appropriate for 
QCD. One main result of the theory concerns the rate of critical slowing down near 
the transition. Just as the correlation length in space diverges at the critical point, 
so does the correlation time for dynamics. Its scaling property is conventionally 
written in terms of a critical exponent z, such that the correlation time scales as 
^^, where ^ is the correlation length. For the model in question, we show below 
that the exponent z governing critical slowing down is predicted to be djl = 3/2. 

Let us now consider how model G of [16] must be modified to deal with the 
four component order parameter of interest. As we discussed above, it is important 
to find the conserved quantities. First, there is the energy. However, it is shown 
in [17] that if a, the specific heat exponent, is negative then the dynamics of the 
order parameter is not affected by the presence of a conserved energy. However, 
the six generators of 0(4) rotations are associated with six conserved quantities 
which can be written as an antisymmetric tensor 



Jap= d Xja(3 ■ (6.1) 

At the level of the unrenormalized Lagrangian, ja can be expressed in terms of the 
order parameter as 

d 

JaP = (^ap-fS(pa-g^(pp , (6.2) 

but this relation need not hold upon renormalization. (Note that since we are 
interested in the critical phenomena, by renormalization we mean renormalization 
towards the infrared.) In model G, the symmetry is 0(3), so there are 3 conserved 
quantities instead of 6. The 3 conserved quantities are the magnetization, while 
the nonconserved order parameter is the staggered magnetization of the antiferro- 
magnet. The equations of motion for j and contain two types of terms. There 
are dissipative terms which damp the system toward the equilibrium configuration 
and so called mode-mode coupling terms. The latter refiect the Poisson bracket 
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relations among the fields. In our case, these are 



0a, Jp'y 



(6.3) 



and 



Jap-, J^S — ^p-fJaS — Sa-yJpS ~ SpsJa-y + ^aSJp-f ■ 



The equations of motion are 



(6.4) 



dt 



J.SF 



0a, J/37 



5F 

Sjp^ 



+ e 



(6.5) 



and 



djap 



7V^ 



SF 



Jap, 07 



SF 



dt Sjap 

where the free energy F is given by 
F = j d^x[- d'rdi<Pa + y 0"0a + j{ 



+ 9 



JaP^ J'jS 



SF 

J Sj-yS 



+ C 



1 

2X^ 



'"X? + T^Japr^ 



H 



a 



(6.6) 



■ (6.7) 



9 and C, are Langevin noise terms. The difference between the dissipative terms 
in the two equations of motion reflects the fact that j is conserved, and hence 
a spatially constant j can not dissipate. Note that j appears in the free energy 
only as p. This is because any higher order terms like j'^cjP' or those involving 
Vj are irrelevant. Because of the form of F, it turns out that the term in the 
equation of motion (6.6) involving (6.4) is zero. The equations of motion given 
above determine the universal dynamics of the long wavelength modes of interest 
near the critical point. 

With equations of motion in hand, Halperin et al. [16] go on to formulate 
dynamic renormalization group transformations. These are obtained by starting 
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with the theory defined with an ultraviolet cut-off A, integrating out modes in the 
momentum shell between b^^A and A and then rescaling according to 



X ^ X = X 



_0/ = b^d-2+r,)/2^ 

(6. 

■lie- 

3^3 = b^3 
t ^t' = h-H . 



After one such transformation, one obtains a new free energy and new equations 
of motion. If one works in the e expansion, these are related to the former F and 
equations of motion simply by a transformation of the parameters /i. A, F, 7, g, 
and X- Repeated application of the renormalization group transformation therefore 
leads to recursion relations for these parameters. The next step is to find the fixed 
point of the transformation. Just as in the static case the fixed point condition 
fixes 77, here it fixes r^, c, and z. Fortunately, to obtain c and z we only need 
the particularly simple recursion relations for x ^-i^d g. Because j only appears 
quadratically in the free energy, the recursion relation for x is 

and this is valid to all orders in e. The terms in the equations of motion proportional 
to g are consequences of the 0(4) symmetry, and as a result the Ward identities 
enforce the recursion relation 

g -^ b'-'^+'g (6.10) 

to all orders in e. From (6.9) and (6.10) it is clear that at a fixed point one must 
have c = d/2 and, as advertised earlier, z = d/2. 
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7. Implications on the Lattice and in the Real World 

Finite temperature lattice QCD simulations are ideally suited to testing many 
of the predictions made in this paper. The static correlation functions of the three 
dimensional theory are natural objects to consider in finite temperature (Euclidean) 
simulations. Also, it is much easier to vary parameters like the temperature and 
the bare quark mass in a lattice simulation than in a real experiment. Hence, it 
should be possible to measure the static critical exponents of section 2, the equation 
of state and the scaling behaviour of the pion and sigma masses of section 3, the 
behaviour of the p and Ai correlation functions of section 4, and the static tricritical 
exponents of section 5. 

Present simulations [14,18] provide strong evidence that the QCD phase tran- 
sition is second order, and that the order parameter is M. [1]. However, there are 
several reasons why present simulations can not yet be used to test the more de- 
tailed predictions of this paper. The fundamental reason is that in all simulations 
to date, the bare quark mass has been so large that correlation lengths do not get 
very long at Tg. For example, in the work of Bernard et al. [18], the correlation 
length in the pion channel at Tc is only about 2.5 lattice lengths. Hence, in order 
to study the behaviour nearer to the critical point and to measure universal prop- 
erties, we must wait for simulations with smaller quark masses. The fundamental 
problem is that long correlation lengths are necessarily accompanied by numerical 
critical slowing down, and this makes simulations challenging. 

Another hurdle to be overcome before lattice simulations can measure the 
critical properties of the QCD phase transition is that any lattice implementation 
of fermion fields only exhibits the full chiral symmetry in the continuum limit. 
If Wilson fermions are used, there is no chiral symmetry at all on the lattice. If 
Kogut-Susskind fermions are used, four flavors of fermions are required, but there is 
a continuous U{1) x f/(l) chiral symmetry on the lattice. This should give a phase 
transition in the universality class of the N = 2 magnet, and has been discussed 
in [19]. In order to study two flavors of fermions, one takes the square root of the 
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fermion determinant in the lattice action. It is not at all clear what this does to 
the lattice chiral symmetries. Hence, for Kogut-Susskind as for Wilson fermions, 
before we are able to test our predictions for the critical phenomena we need finer 
lattices so that extrapolation to the continuum limit can be done. A necessary 
condition for doing this is that the results for Kogut-Susskind and Wilson fermions 
agree when so extrapolated. This condition has not yet been met [20]. 

Hence, while we cannot test our detailed predictions against current lattice 
simulations, we are confident that in the future with finer lattices and longer cor- 
relation lengths, simulations will be able to measure critical exponents, correlation 
functions, and the equation of state, and verify our results. 

Let us now turn to real experiments, as opposed to those on the lattice. Here, 
we are not free to dial the bare quark mass. We will see that this is unfortunate, 
particularly in our discussion of heavy ion collisions. First, however, let us dispose 
of two other possible arenas for testing our results. The QCD phase transition 
certainly occurred in the early universe. Indeed, much work has been done on 
possible observable effects of this transition, if it is first order. Unfortunately, for 
physical values of the quark masses, we have seen that the transition is second order. 
We can think of no observable consequences of a second order QCD phase transition 
in the early universe. In [1], it was noted that since certain antiferromagnets 
including dysprosium have order parameters in the A^ = 4 universality class [21], 
experiments on the phase transition in these materials would help us understand 
the QCD transition. Alas, in these magnets there is a quartic operator like (tt^ + 
•^i)!^! + ^^) which is allowed by the microscopic hamiltonian of the magnets [22], 
but not by that of QCD. This operator makes the symmetric Heisenberg fixed point 
infrared unstable, and either makes the phase transition first order or makes it 
second order but governed by an anisotropic fixed point. Hence, neither cosmology 
nor dysprosium are suitable arenas for learning about the QCD phase transition. 

Now, we turn to relativistic heavy ion collisions. In Bjorken's [23] picture 
of such a collision, a volume of hot plasma forms and quickly reaches thermal 
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equilibrium. In the center of mass frame, the incident nuclei are both Lorentz 
contracted into pancake shapes. They pass through each other, and leave behind 
a region of hot vacuum. In translation, this means that the baryon number of the 
incident nuclei ends up in that part of the plasma heading approximately down the 
beam pipes, and the central rapidity region consists of plasma with approximately 
zero baryon number. In the remainder of this paper, we will attempt to use what 
we have learned about the critical phenomena associated with the chiral phase 
transition to study the behaviour of the plasma in the central rapidity region as it 
expands and cools through T = Tc and eventually hadronizes and becomes pions 
which fly off and are detected. 

The defining characteristic of a second order phase transition is the divergence 
of correlation lengths. How could this feature be observed here? Large volumes 
of space with the order parameter correlated and pointing in a direction different 
than the true vacuum {i.e. sigma) direction will become regions in which the order 
parameter oscillates coherently about the sigma direction. After hadronization, 
correlated volumes will turn into regions of space where the ratio of the number 
of charged pions to neutral pions has some fixed value. Since, in the standard 
scenario of Bjorken [23], different positions in the plasma along the beam direc- 
tion become different longitudinal momenta (actually, different rapidities) as the 
plasma expands, one would hope that a signal of a second order phase transition 
would be fluctuations in the ratio of charged to neutral pions as a function of ra- 
pidity. Coherent evolution of a classical pion fleld has been considered before by 
Anselm and Ryskin [24] and by Blaizot and Krzywicki [25]. However, these authors 
considered neither an equilibrium second order phase transition as we do in this 
section, nor a quench as we will in the next section. 

A prerequisite for the fluctuations discussed above to be observable is that 
the correlation length must get long compared, say, to Tc. To determine whether 
this does indeed happen, we must leave our universality safety net behind, since 
neither Tc nor the magnitudes of correlation lengths are universal. From Figure 2a, 
it is immediately obvious that we have a problem. The longest correlation length 
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is in the pion channel, and the pion mass is increasing with temperature. This 
suggests that the pion mass at T^ is larger than mTy{T = 0) = 135 MeV. This 
is consistent with the fact that in many models the pion mass increases from its 
zero temperature value as the temperature is increased from zero. This result has 
been obtained in chiral perturbation theory using the nonlinear sigma model [26] 
and also for the linear sigma model [27]. Hence, it seems clear that the longest 
correlation length at T^ will be shorter than (135 MeV)^^. This is to be compared 
to Tc itself, which for the case of two massless quarks is around 140 MeV [18]. 
Hence, even though the quark masses are indeed small (~ 10 MeV), the magnetic 
field H proportional to rriq is large enough to prevent any correlation lengths from 
reaching interesting values. 

There is an appropriate quantitative criterion to determine whether a near 
equilibrium second order phase transition leads to dramatic effects. One compares 
the energy in a correlation volume just below Tg with the zero temperature pion 
mass to determine whether or not the correlated volume can become a large number 
of pions. Using current lattice simulations [18], we can make a crude attempt at 
this comparison. The sum of the energy and pressure in a correlation volume is 
about 1/4 the zero temperature p mass. Taken literally, this means that each 
correlation volume becomes only one or two pions in the detector. As we have 
mentioned, current simulations are subject to many caveats, and so this estimate 
should not be taken literally. Nevertheless, it seems clear that the physical value 
of ruq is large enough that in an equilibrium phase transition a correlation volume 
at Tc does not evolve into a large number of correlated zero temperature pions. 
This is not encouraging. 

We pause here for an aside. The reader may be wondering why, when the 
seemingly small equal quark mass rriu = rrid = rriq has such deleterious effects, we 
have completely neglected the difference between the up and down quark masses. 
It was noted in [1] that unequal quark masses allow terms of the type 

AF oc (5m)2(a2 _ Trf + TT^ + 7r|) . (7.1) 

26 



If one is close enough to the critical point that this term matters, one will dis- 
cover an anisotropic fixed point rather than the symmetric Heisenberg fixed point. 
However, while the effect of a common quark mass, namely the mass of the pion, 
is comparable to T^., the effect of (7.1) is much smaller. For example, the QCD 
contribution to the difference in mass between the charged and neutral kaons is 
about 5 MeV [28]. Therefore, we need not worry about (7.1) in real experiments 
since it is much less important than the effect of the "magnetic field" proportional 
to the common up and down quark mass. Of course, (7.1) could be introduced and 
studied on the lattice. 

It seems clear that if the standard scenario for heavy ion collisions in which the 
QCD plasma cools through Tc while staying close to thermal equilibrium is cor- 
rect, then no correlation lengths will get long enough for there to be any dramatic 
observable effects of the phase transition. The chiral "transition" , like the confine- 
ment/deconfinement "transition," will be a smooth crossover. If we were able to 
dial down the quark masses and hence the pion mass, phenomena associated with 
a second order transition would become more prominent. Alas, in the real world, 
unlike on a lattice, we have no such freedom. 

8. Non-equilibrium Phenomena in Relativistic Heavy Ion Collisions 

We are not done yet. The gloomy paragraph with which we ended the pre- 
ceding section began with a conditional sentence. In this section we will consider 
the observable effects in heavy ion collisions if the plasma does not stay close to 
thermal equilibrium through the transition. There are tantalizing hints in cosmic 
ray physics that point in this direction. Among the zoo of high energy cosmic 
ray events known are a particularly peculiar class of events called Centauros [29]. 
These are events with total energy of order 1000 TeV in which many (of order 
100) charged hadrons each with energies of a several TeV and very few photons 
or electrons are seen in a cosmic ray induced shower. In the sample of events in 
there were 5 Centauros, representing about 1% of the events seen with energy 
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of the appropriate order of magnitude. Centauros are peculiar because so many 
charged pions are observed without any of the gammas that would indicate the 
decay of neutral pions. This apparent violation of isospin invariance is puzzling, 
unless one thinks of it in the language of a second order phase transition in which 
these events can be interpreted as the creation of a volume of QCD plasma in which 
the field has fluctuated throughout most of the plasma in some direction in the 
TTi — 7r2 plane. This implies correlation throughout most of a volume of plasma 
large enough that it becomes about 100 zero temperature pions. We convinced 
ourselves above that this could not happen if the plasma remains close to thermal 
equilibrium. Hence, these Centauro events provide a tantalizing hint that it might 
be wise to consider the effects of going from the symmetric phase above Tc to 
the ordered phase rapidly without maintaining thermal equilibrium. This process, 
called quenching, has been much studied in condensed matter physics. 

We noted in the previous section that if the plasma stays close to equilibrium, 
one finds below Tc that the energy in a correlation volume is small compared to the 
zero temperature pion mass. However, since the energy density at temperatures 
well above Tc is much higher than that below Tc, it is reasonable to hope that a 
quench, in which the energy density does not decrease in a quasi-equilibrium fashion 
through Tc, has a better chance of producing correlated volumes of plasma which 
evolve into many zero temperature pions. In this section we will begin to analyze 
the observable effects if relativistic heavy ion collisions proceed via quenching. In 
a real relativistic heavy ion collision, the phase transition will probably occur by 
something in between a slow equilibrium process and a quench. Only experiments 
can determine which description is more appropriate. 

Let us begin by describing more carefully what a quench is, in the context 
of a Heisenberg magnet in greater than two dimensions with no applied magnetic 
field. The system starts at equilibrium at a temperature well above Tc, fluctuating 
among an ensemble of conflgurations with short correlation lengths. One then 
imagines turning the temperature instantaneously to zero. This means that the 
equilibrium conflguration is now an ordered state with the fleld aligned throughout 
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space. However, this is not the configuration in which the system finds itself. 
It is in one configuration from the ensemble appropriate to a high temperature. 
This configuration then evolves according to the zero temperature equations of 
motion {i.e. microcanonically, with no thermal fiuctuations.) In a condensed 
matter system, the appropriate equations of motion are those of Section 6, and in 
particular they are not Lorentz invariant. What is found [30] is that the size of 
correlated domains grows with time in such a way that after a brief initial period 
the correlation function has the simple scaling form 

C(r,t) = (0(r,t)0(O,t)) = (?(r/L(t)) . (8.1) 

The characteristic domain size L increases with time according to L{t) ~ t^, where 
the exponent p depends only on d, N, and the dynamic universality class. It is im- 
portant to note that the scaling behaviour (8.1) is obtained regardless of the initial 
configuration. Hence, it is not actually necessary that the initial configuration in a 
simulation (or in a heavy ion collision) be selected from a high temperature ther- 
mal ensemble. Any disordered initial configuration evolves to the critical behaviour 
(8.1) . This phenomenon is called self-organized criticality. It is also crucial for 
us that the domain size is not related to an equilibrium correlation length, and in 
particular that in an infinite system it grows without bound. 

Let us now consider how the physics of quenching may be applied in relativistic 
heavy ion collisions. If the collision is energetic enough to create a region of plasma 
well above T = T^, the field will indeed be fiuctuating among an ensemble of 
disordered configurations. At the end of the process, one certainly has zero tem- 
perature. The question is what happens to the field in between. One idealized 
possibility which has been considered by many before us and which we considered 
in the previous section is that the system stays arbitrarily close to thermal equi- 
librium. Another idealization is that thermal fiuctuation ceases instanteously. In 
a real collision, the plasma is cooling, and so cannot be exactly in thermal equilib- 
rium. If it cools fast enough, the configuration of the field will "lag," and as in 
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a quench the system will find itself in a configuration that is more disordered than 
the equilibrium configuration appropriate for the current temperature. However, it 
is clear that a real collision will not be an ideal quench, as thermal fiuctuation will 
not cease instantaneously. Thus, a real collision fits in neither idealized category, 
but is somewhere in between. To the second idealization, the quench, we now turn. 

There are several differences between a quench in the system of interest to us 
and in the condensed matter system we described above. First, at zero temperature 
we must have a Lorentz invariant field theory, with different equations of motion 
than those used in the condensed matter system. While the equations of motion 
of section 6 describe the Lorentz non-invariant dynamics of the order parameter 
near T = Tg, they are inappropriate for the T = dynamics of a quench. We 
propose to use the zero temperature linear sigma model with only pion and sigma 
fields. Since we will be considering energies well below the rho and nucleon masses, 
we need not include these degrees of freedom. The sigma could also be left out, 
if it were not for the fact that the appropriate initial conditions are a disordered 
state in which the pion and sigma fields are equivalent up to the effect of the bare 
quark masses. The second difference is that unlike in the magnet, the plasma in 
a relativistic heavy ion collision is expanding. This means that the description 
in terms of configurations of the field will not be appropriate forever. At some 
time, the energy density drops low enough that one has individual pions fiying off 
towards the detectors. The third difference is that unlike in the condensed matter 
systems considered in [30] , we must include the effects of the bare quark masses. 
As we saw before, these correspond to a significant magnetic field. We propose 
that quenching of a 4 component Heisenberg model with the three modifications 
we have mentioned be considered an idealized model for a heavy ion collision. 

It is fortunate that the scenario we have just outlined, with the exception of the 
significant magnetic field, is exactly the scenario considered by Turok and Spergel 
[31] as a cosmological model for large scale structure formation in the early universe. 
They study the evolution of an 0{N) sigma model in an expanding universe. They 
find an exact scaling solution for the non-linear sigma model in the large A^ limit, 
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and do numerical simulations for the linear sigma model for A^ = 4 and A^ = 10. 
The main reason we can apply their results while we cannot use those of [30] is 
that Turok and Spergel use the Lorentz invariant equations of motion appropriate 
for our problem. They find that the size of correlated domains, L{t), grows at the 
speed of light! We are currently [32] extending their simulations to include the 
effects of a magnetic field, and to vary the expansion rate. The magnetic field will 
qualitatively change the scenario. Instead of having correlated domains with the 
field pointing in arbitrary directions on the 3-sphere, at late times the field 
will be oscillating about the sigma direction. However, there will still be domains 
in which the oscillations are in different directions. If in a heavy ion collision the 
size of these domains grows with the speed of light as it does in our simulations 
[32], the phenomenological consequences for heavy ion collisions are dramatic. Of 
course the description we are using will only be valid for a short time. (Bjorken [23] 
estimates that a hydrodynamic description will be valid for about 10 fm/c.) Since 
the plasma will be expanding slower than the speed of light, even in this short 
time domains which expand at the speed of light will grow to encompass large 
fractions of the total volume. We therefore propose that if heavy ion collisions can 
be modelled as a quench, this will be detected by observing clusters of pions in 
which all the pions in a region of rapidity are correlated in internal space. In some 
clusters, there will be only charged pions; in others, only neutral ones; and in all, 
charged and neutral pions will occur in some fixed ratio. 

We can estimate the probability distribution of the ratio R of the number of 
neutral pions to the total number of pions in a correlated region. Let us assume 
that the field in the region is initially equally likely to be pointing in any direction 
on the 3-sphere. This assumption may not be strictly true because the magnetic 
field selects a preferred sigma direction even at high temperatures. However, we 
make the assumption in order to get a simple analytical result. Where starts 
will determine in which direction it ends up oscillating about the sigma direction. 
We define angles on the 3-sphere according to 
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(a, TTs, TTi, 7r2) = (cos^, sin6'cos0, sin^sin^cosr], sin6'sin0sinr7) . (8.2) 

Then the ratio R is given by 



sin cos^ 2 

n^o + n^+Tj- sin^ ^(cos^ cp + sin^ 



it = = TT-, 7. 9-— = COS . (8.3) 



Under the assumption that all initial values on the 3-sphere are equally probable, 
the probability distribution V{R) is determined by 

R2 2-K -K arccos(\/^) 

I V{R)dR = \ I dr] I ddsiii^d I rf0sin0 (8.4) 

■Ri 00 arccos(/Rl) 



V{R) = \r'^''^ . (8.5) 



and turns out to be simply 



Equivalently, the probability that i? < -Ri is given by yRi. If heavy ion collisions 
are described by a quench, there should be large regions of the collision volume 
containing clusters of pions in which R is constant, and the values R takes in 
different such regions should be distributed according to (8.5) . 

As one application of (8.5) , we note that the probability that the neutral pion 
fraction R is less than .01 is 0.1! This is a graphic illustration of how different (8.5) 
is from what one would expect if individual pions were independently randomly 
distributed in isospin space. It also makes Centauro events in which less than 1% of 
the outgoing particles from a heavy ion collision are neutral pions seem much less 
surprising than they first appeared. The analysis of the Centauro data is difficult 
for several reasons. Most important of these is the limitation imposed by small 
statistics. Also, if a Centauro event occurs too high above the detector, so many 
secondary photons will be produced that the event will not be recognized as a 
Centauro. Third, in a Centauro event all of the particles from the collision strike a 
small region of the detector and it is impossible to isolate the central rapidity region. 
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This combined with the fact that the detectors do not distinguish between charged 
pions and charged baryons has several unfortunate consequences. It means that 
events in which there are two or more correlated clusters of pions are not detected 
as Centauros. Only those with a single cluster are so identified. Also, the opposite 
type of event in which R is close to 1, will not have a dramatic signature since there 
will always be charged baryons present from the two initial nuclei. For all these 
reasons, we feel it is impossible to extract a meaningful probability distribution 
'P{R) from the cosmic ray data. 

When relativistic heavy ion collisions occur at high enough energies in a labora- 
tory colliding beam facility, all of the difficulties of the cosmic ray experiments will 
be rapidly overcome. That will be the time to look for correlated clusters of pions, 
and to look for a distribution like (8.5) , and hence to determine whether these 
collisions proceed by a process close to the idealized quench we have considered 
here. 

9. Conclusions 

The future of the study of the QCD phase transition on the lattice looks promis- 
ing. As simulations improve, they will begin to investigate the plethora of static 
critical phenomena we discussed in the first sections of this paper. Critical expo- 
nents, the equation of state, the critical behaviour of the pion and sigma suscepti- 
bilities, p and Ai correlation functions, and tricritical exponents are all out there 
waiting to be measured. 

Because heavy ion collisions are dynamical processes, there are more possible 
scenarios and the situation is not clear cut. We have considered two idealized 
models. If the cooling plasma stays arbitrarily close to thermal equilibrium, our 
conclusions are disappointing. Because the pion is so heavy compared to Tg, corre- 
lation lengths will not become particularly long. On the other hand, in the other 
idealized model we considered, in which thermal fiuctuations are rapidly quenched 
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and the cooling plasma is in a sense maximally out of thermal equilibrium, dra- 
matic phenomena are possible. Correlated volumes will form, and their size will 
not be determined by any equilibrium correlation length. Indeed they may grow at 
the speed of light until hadronization occurs. This will have the consequence that 
clusters of large numbers of pions will be detected in which the ratio of neutral to 
charged pions will be constant. This ratio will be different in different clusters, and 
will follow a probability distribution which is skewed towards having few neutral 
pions. We eagerly await the verdict of experiment as to which scenario is more 
appropriate. 
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APPENDIX 

In this appendix, we reproduce the equation of state to order e^ [9] , and give 
various other results used in section 3. 

The equation of state is found by doing an e expansion of the relation 

{a) = 0, (A.l) 

where a = a — M. The resulting expansion can be expressed in terms of the 
variables y = H/M and x = t/M^'" ^ where 



13 AT + 8 (AT + 8^3 



^ = 2 + -TT^e + 4 ^^V.. ^ r. ^ ^' + 0{e') (A.2) 



and 



A^2 + UN + 60 2 

2(A^ + 8^2 



3\ 



(5 = 3 + 6+ '.^^ ^,; 6^ + 0(e^). (A.3) 



We are interested in e = 1 and A^ = 4. We choose to measure fields in units such 
that y = 1 at t = 0, and a; = — 1 at the coexistence curve {H = 0, t < 0). The 
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equation of state is 

y = fix) . (A.4) 

The function f{x) is given by 
f{x) = 1 + a; 

+ 2(ivT8) (' + 2(ivT8) [^ -' + «'" 2 - '"° ^ + '^ - ^""'^ +')] ) 

X [3(a; + 3) ln(a; + 3) + (A^ - l)(a; + 1) ln(a; + 1) + 6a;ln2 - 9{x + 1) ln3) 

+ 36 [In^ (x + 3)-{x + l)\n^3 + x In^ 2] 

27 

- 54 In 2 [ln(a; + 3) + a; In 2 - (x + 1) In 3] + 3(A^ - 1) ( In — ) {x + 1) ln(a; + 1) 

212 + 17N - 4iV2 
+ — [{x + 3) ln(x + 3) + 2a; In 2 - 3(a; + 1) In 3] 

+ (iV - l)(a; + 1) ln(a; + 1) ln(a; + 3) {N - l){x + 1) ln2(a; + 1) 

+ ]y^(19^ + 92)(a; + 1) ln(x + 1) - 2{N - 1) [{x + 6)/i(p) - Q{x + l)/i(3/4)] 

- 6{N - 1) [hip) -{x + l)/2(3/4)] + 4(A^ - 1) [hip) -{x + l)/3(3/4)] } + 0{e^) 

(A.5) 



where 



and 



4(x + 1) 



oo 

dulnu 



J u{l-u)^ 'J u{l 



u 



(A.7) 



hip) = p-r- 
dp 

hip) ^ hip) + 2hip) . 
The behaviour of h (p) near the coexistence curve is given by 

-^i(p) ~ -;-(ln^4p + 21np) , p ^ oo . (Ai 

4p ^ 
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The leading terms of f{x) for large x are 



where 



^(")~(^ + I^W'''^^''^^"'' "^" 



^ ^ I ^ + ^ , I (A^ + 2)(iV^ + 22iV + 52) 2 , ^.3 



(A.9) 



(A.IO) 



Of course, this result for 7 and the results for (3 and 6 in (A2) and (A3) are 
consistent with the scaling relations (2.12) . We saw in section 3 that the behaviour 
(A.9) determines the sigma and pion masses for iJ = and t > 0. 

In section 3, in order to determine m^ near the coexistence curve we needed 
the behaviour of f'{x) for x ^ —1. In this region, f{x) is given by [11] 



fix)r^{x + l){l + e 



^-1 w .N 3(l + 31n(2/3)) 
■ln(a; + l) + ^ ^ ' '' 



2{N + 



2{N + 8) 



f^-i^'f^-i'm^f. + i; 



+ 



8(A^ + 8)2 
A^-1 



4(iV + 



iV + 27+181n2-91n3 



60 



iV + SJ 



One next inverts (A. 11) to obtain x + 1 in terms oi y = f{x), 



(a; + l)=Cit/ + C2y^-^/2 + 0(l/2-0(^ 



ln(a; + 1 

(A.ll) 

(A.12) 



O(e') 



differentiates the result with respect to x, and obtains 

1 



^^ -. f'(x) 



d + C2y ^/2 



for X 



(A.13) 



with 



ci 



N 



2{N- 



(iV + 8)ln2-91n3 + 



25iV2 + 142A^ + 76 
9(iVT8) , 



(A.14) 
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and 



C2 = C2 ( 1 - - 



N-1 e , 

l-:r 1 



-^—2 [9{N + 8) In 3 + 22iV + 116] U + 0{e'^ 



N + 8[ 2 V {N 

(A.15) 
As H is lowered to zero at fixed t < 0, m^ at first tends toward a constant, and 
then goes to zero according to m^ oc H^'^ when the C2 term takes over from the 
ci term. For A^ = 4, this occurs when 

H /C9\2/e /l \2/e 

W^icJ -(gCl + O-See]) ~0.4fore=l. (A. 16) 

Of course, the numerical value for e = 1 should not be taken too seriously. The 
qualitative result is clear nevertheless. 

We end this appendix by describing how the figures in section 3 were obtained. 
When evaluated at e = 1, the expression (A. 5) for f{x) has several problems. First, 
at large x it does not grow as x^ . Rather, it increases like x In x. Also, for x -^ —1, 
f{x) given by (A. 5) does not satisfy (A. 12) . In fact, for x ^ —0.95, f{x) < 
which is unphysical. Both of these problems arise because we are setting e = 1 in 
an expansion of f(x) to order e^ which is valid for e ^ 0. In order to illustrate 
the correct qualitative behaviour, we constructed a function f{x) which smoothly 
interpolates between (A. 9) at large x, (A. 12) near x = —1, and (A. 5) in between. 
Using this function f{x), we obtained Figure 1 by solving H/M = f{t/M^'") 
for M at various values of t and H. We then calculated the results for m^ and 
m^ shown in Figures 2 and 3 using (3.6) and (3.8) . Because of the limitations 
imposed by working at e = 1, all three figures should be viewed as illustrations of 
qualitative behaviour. 
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FIGURE CAPTIONS 

1) The order parameter M as a function of reduced temperature t = {T — Tc)/Tc 
for magnetic fields if = 0, 0.002, 0.005, and 0.02. Like t, both M and H are 
dimensionless. They are obtained from their dimensionful counterparts by 
dividing by non-universal dimensionful constants defined in such a way that 
for t < and H = the order parameter is given by M = (— t)^, and for 
t = it satisfies M = H^/^. 

2) ml (Figure 2a) and m^ (Figure 2b) as functions oft for H = 0, 0.002, 0.005, 
and 0.02. Since M and H are in dimensionless scaled units, so are m^ and 
m^. For t = 0, m^ = H^^^^^/^ and m^ = 5rn^. For H = and t > (and for 
large enough t for any H) m^ = m^ ~ f^. For t < and /J — > 0, m^ ~ H 
and the sigma mass decreases to zero as shown in Figure 3. 

3) m^ as a function of H for t = —0.2. For large H it behaves as if it will be 
non-zero at H = 0, but in fact for if — * it decreases like m^ ~ H^ where 
to lowest order in e, p = e/2 = 1/2. 
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